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Abstract 

The Damage Spreading (DS) method allows the investigation of the effect caused by tiny perturbations, in 
the initial conditions of physical systems, on their final stationary or equilibrium states. The damage (D(t)) 
is determined during the dynamic evolution of a physical system and measures the time dependence of the 
difference between a reference (unperturbed) configuration and an initially perturbed one. In this paper we first 
give a brief overview of Monte Carlo simulation results obtained by applying the DS method. Different model 
systems under study often exhibit a transition between a state where the damage becomes healed (the frozen 
phase) and a regime where the damage spreads arriving at a finite (stationary) value (the damaged phase), when 
a control parameter is finely tuned. These kinds of transitions are actually true irreversible phase transitions 
themselves, and the issue of their universality class is also discussed. 

Subsequently, the attention is focused on the propagation of damage in magnetic systems placed in confined 
geometries. The influence of interfaces between magnetic domains of different orientation on the spreading of the 
perturbation is also discussed, showing that the presence of interfaces enhances the propagation of the damage. 
Furthermore, the critical transition between propagation and nonpropagation of the damage is discussed. The re- 
sults analyzed indicate that, in some cases as in the Abraham's Model and in the standard Ising magnet (Glauber 
dynamics), there is clear evidence showing that the DS transition and the critical transition of the physical system 
(in these cases the wetting and the ferromagnet-paramagnet transitions, respectively) occur at different critical 
points. However, in the case of the corner geometry, the critical points of both transitions— damage spreading 
and cornerfilling— coincide within error bars. It is found that, at criticality, the damage obeys a power-law 
behavior of the form D(t) oc t v , where r\ is the damage spreading critical exponent. The evaluation of critical 
exponents allows the identification of three propagation regimes: i) inside the magnetic domains the propagation 
is slow [r\ — 0.40(2)), ii) the fast propagation is observed along the interface between domains (rj = 0.90(2)), iii) 
the alternating propagation across interfaces and inside domains is consistent with an exponent lying between 
the previous cases, namely 77 = 0.47(1). In all cases, the determined critical exponents suggest that the DS 
transition does not belong to the universality class of Directed Percolation, unlike many other systems exhibit- 
ing irreversible phase transitions. This result reflects the dramatic influence of interfaces on the propagation of 
perturbations in magnetic systems. 

Keywords: Propagation of perturbations; Phase transitions and critical phenomena; Magnetic systems in 
confined geometries. 
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1 Introduction. 



One of the most interesting challenges in the theory of dynamic systems is the understanding of the dependence 
of the time evolution of physical observables on the initial conditions, because very often a small perturbation in 
the initial parameters could completely change their behavior [1]. Within this context, it is interesting to study 
the time evolution of such perturbations in order to investigate under which conditions a small initial perturbation 
may grow up indefinitely or, eventually, it may vanish and become healed. 

In order to understand this behavior, Kauffman introduced the concept of Damage Spreading (DS) [2]. In order 
to implement the DS method in computational simulations [3,4], two configurations or samples S and S' , of a 
certain stochastic model, are allowed to evolve simultaneously. Initially, both samples differ only in the state of a 
small number of sites. Then, the difference between S and S' can be considered as a small initial perturbation or 
damage. 

The time evolution of the perturbation can be followed by evaluating the total damage or "Hamming distance" 
defined as 

i i 

where Di(t) is the damage of the site labeled with the index i at time t, £s 4 (t),s<(t) is the delta function and the 
summations in equation ([I]) run over the total number of sites of the system N. 

By starting from a vanishing small perturbation D(t = 0) — * 0, one may expect at least two main scenarios, 
namely: a) D(t — ► oo) — > and the perturbation is irrelevant because the damage heals; or b) D(t — ► oo) assumes 
some non-zero value, and the damage spreads. These two situations indicate a weak or a strong influence of the 
initial conditions, respectively [5,6]. Of course, an intermediate or marginal case where the number of damaged sites 
approaches a finite positive number, which does not increase proportionally to N, may also occur. In this case, after 
normalization by N one has that the equation (TJ gives D(t — > oo) = 0, but the damage has not healed out. When 
the system arrives either at the state of zero damage or complete damage (D = or D = 1, respectively), it will 
remain in this situation indefinitely. For this reason, a transition between these two states is irreversible and could 
be related to Directed Percolation processes [7], for which an irreversible critical transition occurs from an active 
to an inactive state (absorptive state). In those models belonging to the universality class of directed percolation, 
when the system is trapped in an absorptive state, it is then impossible to recover the activity by changing the 
control parameter. 

The first studies of DS in physical systems were applied to the Ising model, spin glasses and the Kauffman cellular 
automata, and they appeared in the mid-eighties [8-11]. Subsequently, this technique has also been applied to the 
study of several different models such as Ising models [9-59], spin glasses [11,60-70], Potts models with q-states 
[71-77], the Heisenberg model [78,79], the XY model [80,81], a discrete ferromagnet [82], two-dimensional trivalent 
cellular structures [83,84], biological evolution [85,86], and cellular automata [6,87-93]. Also, non-equilibrium 
systems [2,28,92,94-99], SOS models [100-102], opinion dynamics [103,104] and small world networks [105,106], 
have been characterized by means of the DS method. 

In view of the great interest attracted by this field of research of interdisciplinary application, the aim of this 
work is to present a brief overview of the state of the art in studies of the Damage Spreading transition, focusing our 
attention on recent results obtained for magnetic systems in confined geometries, which are helpful in understanding 
the propagation of perturbations in nano- and micromaterials. 

The manuscript is organized as follows: in Section 2 we briefly describe the archetypical models used in the 
study of DS, namely the Domany-Kinzel cellular automata and the Ising model. In Section 3 we describe the DS 
method in these basic models, and Section 4 is devoted to the discussion on the main characteristics of DS. In 
Section 5, in first place we present a brief discussion of the equilibrium configurations that can be found for Ising 
systems in confined geometries, and therefore we discuss the results obtained for DS in these systems. Finally, our 
conclusions are stated in Section 6. 

2 Definition of basic models. 

Both, the Domany-Kinzel (DK) cellular automata model and the Ising magnet have become archetypical systems 
for the study of DS. So, for the sake of completeness we give brief descriptions of both models. 
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2.1 The Domany-Kinzel Cellular Automata Model. 



The DK model [107, 108] is a family of the (1 + 1) dimensional stochastic cellular automata with two parameters, 
Pi and P2, which simulate the time evolution of interacting active elements in a random medium. The DK model 
consists of a linear chain of N sites (Si) that can take two possible values, usually and 1 (empty and occupied 
sites, respectively). The state of each site i at time t + 1 [Si(t + 1)] depends only upon the state at time t of the 
two nearest neighbors [Si-±(t) and Sj+i(t)], according to the transition probability [P[Si(t + 1) | Si-i(t), Si+i(t)]] 
defined as 

P[1|0,0] = (2) 
P[1|0,1] - P[l|l,0]=pi 

^[1 I 1,1] = P2 , 

where P[0|»,»] = 1 — P[l|»,»] and the parameters p\ and pi represent the probabilities that the site i is occupied 
if exactly one or both of its neighbors are also occupied, respectively. 

Domany and Kinzel demonstrated [107, 108] the existence of two phases, depending on the values of the pa- 
rameters pi and P2, a frozen and an active phase, separated by a critical line. In the active phase, there exists a 
stationary state that it is governed by fluctuations, while in the frozen phase all initial states lead to an absorbing 
state. There is strong numerical evidence that this phase transition belongs to the universality class of Directed 
Percolation (except for its upper terminal point). The exceptional behavior at the upper terminal point of the 
critical line is due to an additional symmetry between active and inactive sites along the line P2 = 1. Here, the DK 
model has two symmetric absorbing states given by the empty and the fully occupied lattices, respectively. 

2.2 The Ising Model. 

The other archetypical system used to study the DS transition is the Ising model [109]. In this case, each site of 
the lattice represents a spin variable. In the ferromagnetic case, the spins have an energetic preference to adopt the 
same direction. The Hamiltonian of this system can be written as 

n = -j- oiOj- H Yl ai ( 3 ) 

<i,j> <i,j> 

where <Tj is the Ising spin variable that can assume two different values <Ji — ±1, the indexes 1 < i, j < N are used 
to label the spins, J > is the coupling constant of the ferromagnet, H is the external magnetic field, and the 
summation runs over all the nearest-neighbor pairs of spins. In the absence of an external magnetic field (H = 0) 
and at low temperature, the system is, for more than one dimension, in the ferromagnetic phase and, on average, 
most spins are pointing in the same direction. In contrast, at high temperature the system maximizes the entropy, 
thermal fluctuations break the order and the system is in the paramagnetic phase. This ferromagnetic-paramagnetic 
critical transition is a second-order phase transition and it occurs at a well-defined critical Temperature (Tc). In 
the two-dimensional case, one has exactly kTc/J = 2/ln(l + \/2) = 2.269...., where k is the Boltzmann constant. 

3 Damage Spreading in the basic models. 
3.1 Damage Spreading in the DK Model. 

One of the most interesting behaviors of DS was first discovered in the (1 + l)-dimensional DK cellular automata 
[107, 108]. In fact, in addition to the two known phases (frozen and active) of the phase diagram of the DK model, 
already discussed in the previous section, Martins et al. [110] found a third phase related to the spreading of 
the damage. This "new" phase, which lies in the active phase, simultaneously exhibits regions where the damage 
spread and heals, and it is extremely sensitivity to the initial conditions. Subsequently, other authors [34,111-113] 
determined the boundary of this phase more precisely. Independently, Mean Field approximations applied to 
different systems [90,93,113,114] confirmed the existence of this "chaotic phase" in the sense that the final state of 
the system is sensitive to the initial conditions. However, it has also been realized that the results obtained [93,114] 
may depend on the dynamic rules used in the implementation of the algorithm. This characteristic of DS will be 
explained in detail in next Section. 
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3.2 Damage Spreading in the Ising Model. 

In the case of the Ising model, the definition given in equation (TJ can be rewritten as 

1 N 

D(t) = — y £\S?(t,T)-S?(t,T)\, (4) 
/ 

where the summation runs over the total number of spins N, and the index I (1 < I < N) is the label that identifies 
the spins of the configurations. S A (t, T) is an equilibrium configuration of the system at temperature T and time 
t, while S B (t,T) is the perturbed configuration that is obtained from the previous one, at t = 0, just by flipping 
few spins [3] . Physically, the definition given by equation (@]) represents the total fraction of spins that are different 
in both configurations. Then, one is interested in investigating under which conditions a small initial perturbation 
will grow up indefinitely or eventually will vanish and become healed. 

Notice that by performing Monte Carlo simulations, the same sequence of random numbers has to be used for 
both copies, in order to assure identical realizations of thermal noise. Furthermore, to study the time evolution of 
the DS, a meaningful definition of the Monte Carlo time step (mcs) is necessary. For this purpose, the standard 
definition is adopted according to that during one mcs all N— spins of the sample are flipped once, on the average. 

4 Main characteristics of the Damage Spreading. 

4.1 On the Dependence of Damage Spreading with the Simulation Algorithm. 

As it was mentioned previously, the DS behavior depends on the dynamic rules used to implement the algorithm, 
in Monte Carlo Simulations. In order to understand this dependence, it is very useful to remind the reader that the 
detailed balance condition [1,4] assures that the system will arrive at an equilibrium state, but it does not establish 
the way of this evolution. In other words, detailed balance does not univocally determine the dynamic rules that 
have to be applied to go from a given configuration to the next. Therefore, this situation opens the possibility on 
choosing different transition probabilities between states, which implies that different dynamic rules can be applied 
to the same physical system. In the case of the Ising Model, the equilibrium state can be generated by different 
dynamic rules, e. g. Heath Bath, Glauber, Kawasaki and Metropolis dynamics, which represent different dynamics 
that allow the system to arrive at an equilibrium state. 

In principle, it was expected that DS would not depend on the intrinsic dynamics of the simulation, and one 
could find regular and chaotic phases that could be identified with the properties of the equilibrium system [10,11]. 
However, subsequent studies showed that different dynamic rules implies the occurrence of different behavior in the 
DS, such as Glauber versus Metropolis [17,115], Q2R [16] or Kawasaki [19]. Furthermore, depending on the type 
of updates, this means in which way the sites of the lattice are chosen (randomly, typewriter, chessboard, etc.), the 
results obtained for DS and the critical temperature of the transition between propagation and nonpropagation of 
the damage are different [20,24,46]. 

Summing up, the observed behavior is not surprising since DS is a dynamic process, and for this reason, it is 
reasonable to expect that it may depend on the dynamic rules applied. As an example, it is useful to consider 
the Ising model in two dimensions and compare the results obtained using both Glauber [10] and Heath Bath [11] 
dynamics. For the Ising model simulated by using the Glauber dynamics, Stanley et al. [10] and Mariz et al. [17] 

found that in the paramagnetic phase (T > Tc, where Tc is the Onsager critical temperature (Tc — 2.269 J)) 

the system is chaotic (D(t — > oo) tends to a finite non-zero value), while in the ferromagnetic phase the damage 
heals (D(t — > oo) — > 0). On the other hand, Derrida et al. [11] studied DS in the Ising Model with the Heath 
Bath dynamics. They found, in contrast to the case of Glauber dynamics, that for both the paramagnetic and 
ferromagnetic phases, damage heals in the limit D(0) — > 0. 

On the other hand, the results obtained by U. Costa [13] and Le Caer [14] in three dimensions and using the 
Glauber dynamics, show a behavior similar to the 2D case, but the critical temperature for the DS transition (Td) 
is lower than Tc. 

4.2 On the Universality Class of the Damage Spreading Transition. 

The universality class of the continuous and irreversible critical transition between propagation and non-propagation 
of damage is still an open question. Grassberger [34] conjectures that the DS transition may belong to the Directed 
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Percolation universality class (DP) [116-118] if its critical point does not coincide with a critical transition of the 
physical system, e.g. the critical temperature of the Ising magnet. In the same paper, Grassberger [34] presents a 
Monte Carlo simulation study of DS in the DK cellular automata [119] as a test for his conjecture. He founds that 
critical exponents of the DS transition coincide, within error bars, with the exponents of the DP universality class, 
both in two and three dimensions. 

An analytical justification for Grassberger's conjecture is given by Mean Field studies reported by Bagnoli [93] 
and the subsequent exact results obtained by Kohring and Schreckenber [114]. They found that, for certain limits 
of the transition probability between states, the dynamics of DS in the DK cellular automata is identical to the 
evolution of the DK itself, and for this reason, the DS transition belongs to the DP universality class. These results 
were later extended to other regions of the phase diagram of the DK Model [91]. 

Numerical simulations of different models also showed that the DS transition is characterized by critical expo- 
nents of the DP universality class, as in the case of the 2D Ising Model with Swendsen-Wang dynamics [48], as well 
as in a deterministic cellular automata with small noise [87]. 

However, there are also other systems where the DS transition has a non-DP behavior, such as the case of the 
Kauffman Model [2]. In general, it is expected that if the order parameter has a Z2 symmetry, the DS transition 
may not belong to the DP class [40]. 

In the Ising model with Glauber dynamics, Stanley et al. [10] determined that the DS transition and the 
paramagnet-ferromagnet transition coincide. For this reason, it is expected that the DS transition may not belong 
to the DP universality class. 

In three dimensions, Costa [13] and Le Caer [14] found that Tp/Tc ~ 0.96 and Td/Tq ~ 0.91, respectively. A 
few years later, Grassberger [33] determined more precisely the critical temperature for the DS transition, and he 
found that Tjj — 0.992(2)Tc in two dimensions and Td = 0.9225(5)Tc in three dimensions. He also determined the 
value of the critical exponent 6 governing the time dependence of the survival probability of damage, showing that 
it coincides, within error bars, with the accepted value corresponding to the DP universality class. 

On the other hand, Vojta [41,43], studying the dynamic stability of the kinetic Ising model with Glauber 
dynamics and using a Mean Field approximation, found that there exists a critical temperature for the Damage 
propagation-nonpropagation transition given by Td — 1.739J ~ 0.826Tc- Also, Vojta generalized these results in 
the presence of an external magnetic field h. He found that h causes an increase of the critical temperature Td and 
stabilizes the non-chaotic phase. The general behavior of Td as a function of h is given by 

in the limit h/J — > 1. (5) 

DS has also been studied in the Ziff-Gulari-Barshad (ZGB) model [120], for the catalytic oxidation of carbon 
monoxide. The ZGB model exhibits a second-order irreversible phase transition between an active state with 
production of CO2 and a poisoned (absortive) regime [120] where the reaction stops irreversibly, which is known 
to belong to the universality class of DP [121]. By performing Monte Carlo simulations of DS in the ZGB model 
in 2D, Albano [95] showed that there exist both chaotic and regular phases, and that the transition between them 
lies within the reactive phase. So, the DS transition is not coincident with the model transition that takes place 
between a poisoned and a reactive phase. However, the reported critical exponent S ~ 0.65(2) is quite different from 
that corresponding to the DP universality class, namely 5 — 0.451 [118]. So, the question about the universality 
class of the DS transition in the ZGB model is still open. 

Another interesting possibility is to study the interplay between reversible (equilibrium) critical transitions of 
confined systems and the irreversible DS transition. In particular, for the case of magnetic systems confined between 
rigid walls, it is interesting to study the relationship between on the one hand, the DS transition and on the other 
hand, the paramagnet-ferromagnet, the wetting and the corner-filling transitions. In the last examples (wetting 
and corner transitions), the presence of external magnetic fields applied to the walls of the system promotes the 
presence of interfaces between magnetic domains of different orientation, and therefore it is interesting to study the 
effect of these interfaces on the propagation of damage. In fact, it have been observed [54-56,59] that the presence 
of interfaces enhances the spatiotemporal spreading of damage. 

5 Damage Spreading in the Confined Ising Model. 

Here, we focused our attention to the study of DS in the Ising Model when the ferromagnet is confined between walls 
that exhert surface magnetic fields. For this purpose, first it is worth studying the main properties of equilibrium 
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configurations, since they are the starting point for the study of DS. Subsequently, the study of DS is actually 
addressed. 



5.1 Equilibrium Configurations of the Ising Magnet in Confined Geometries. 
5.1.1 Strip Geometries and the Wetting Transition. 

The study of the properties of thin— confined— films has attracted growing attention in the last decades not only 
due to the interest in the understanding of the properties of confined and low dimensional materials, but also to 
the existence of many potential applications in the fields of nanoscience and nanotechnology. 

A thin film can be modeled by using the Ising magnet in a confined geometry. So, if one has a film in a stripped 
geometry of size L X M (L <C M) , the Hamiltonian of the Ising magnet given by equation ([3]) can be redefined as 



where is the Ising spin variable corresponding to the site of coordinates J > is the coupling constant of 

the ferromagnet and the first summation of iJH) runs over all nearest-neighbor pairs of spins such as 1 < i < M and 
1 < j < L. The second (third) summation corresponds to the interaction of the spins placed at the surface layer 
j = 1 (j = L) of the film where the surface magnetic field hi (h^) acts. Such fields are measured in units of the 
coupling constant J. Of course, for non-vanishing surface fields one has to assume open boundary conditions (OBC) 
along the M— direction of the film. However, for h% = Hl = both OBC and periodic boundary conditions (PBC) 
may be used. Also, along the L— direction of the strip, OBC are always used. 

Considering PBC connecting the upper and lower surfaces of the film and neglecting surface fields in the Hamil- 
tonian of equation ((6]), one finds standard configurations of the Ising magnet, as is shown in Figure Q] (a). In this 
case, for T — 0.987c one has rather homogeneous configurations essentially showing large monodomain structures. 
On the other hand, by assuming OBC and keeping hi = hi, = 0, a quite distinct behavior is observed, as is shown 
in Figure Q] (b) , which was also obtained at the same temperature as in Figure Q] (a) . In fact, near the bulk critical 
temperature, the system shows a quasi-phase transition from a state that is ordered at scales £ < L (note that £ is 
the standard correlation length) below Tc to a state that is essentially ordered in the direction perpendicular to the 
open boundary, but is disordered in the other direction. In fact, for T < Tc (but close to Tc) the system is broken 
up in a sequence of magnetic domains with spins of opposite sign. The successive domain walls occur essentially at 
random (see e.g. Figure Q] (b)). It should be noticed that this particular kind of configuration is the macroscopic 
manifestation of the mixing neighbor effect undergone by spins at the surfaces of the film with open boundaries. 
For a detailed discussion on this kind of configuration, see also Albano et al. [128]. 

Also, by considering OBC and competing surface magnetic fields hi — —hi, (see equation ([6])) another very 
interesting scenario takes place. In fact, for T < Tc, these competing fields cause the development of a domain 
wall interface along the direction parallel to the surface of the film, as is shown in Figure [21 This situation can be 
described in terms of a wetting transition that takes place at a certain— field-dependent— critical wetting temperature 
T w (h). In fact, for T < T w (h) a small number of rows parallel to one of the surface of the film have an overall 
magnetization pointing to the same direction as the adjacent surfaces field (see Figure [2] (a)). However, the bulk of 
the film has the opposite magnetization (i.e. pointing in the direction of the other competing field). Alternatively, 
one may also consider a symmetric situation that is equivalent to the previous one due to the spin-reversal field- 
reversal symmetry. This non-wet state of the surface takes place at low enough temperatures. As the temperature 
is raised towards T w (h), the number of rows adjacent to the surface that has a magnetization of different orientation 
than the bulk increases, i.e. the domain wall between the coexisting phase of opposite magnetization, which at 
low temperature is tightly bound to one surface (non-wet state), moves farther and farther away from the surface 
toward the bulk of the film. When the interface is located, on average, in the middle of the film, the system reaches 
the wet phase for the first time. Of course, a well-defined wetting transition takes place in the thermodynamic 
limit only. Nevertheless, as shown in Figure EJb), a precursor of this wetting transition can also be observed in 
confined geometries for finite values of L. In confined geometries, this precursor wetting transition is most correctly 
described in terms of a localization-delocalization transition of the interface (for further discussions on the wetting 
transition of the Ising, system see e.g. [122-137]). For T > T w , the interface between domains moves along the 
i-direction, and the system enters the wet regime (delocalized interface), as is shown in Figure [2] (c) . 
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The phase diagram (i.e. the critical curve in the h — T plane, as shown in Figure \3§ has been solved exactly by 
Abraham [122,126], yielding 

cosh(2fy3) = cosh(2if ) - e- 2K sinh(2JC), (7) 

where J > is the coupling constant, h is the surface magnetic field, /3 = 1/kT is the Boltzmann factor, and 
K = J/3. 

5.1.2 The Corner Geometry and the Filling Transition. 

Another interesting confinement scenario for the Ising magnet is the so-called corner geometry sketched in Figure 
IH For this case, the Hamiltonian given by equation ^ has to be modified for a lattice of size L x L, yielding 

H = -J Vi,jPm,n — h cr^i — h y~] cxlj + h a-jj + h a it L, (8) 

<i, i 7,m,n> i j j i 

where <7ij = ±1 is the spin variable, J > is the coupling constant, and h is the magnitude of the surface field. 
The first summation runs over all spins, while the remaining ones hold for spins at the surfaces where the magnetic 
fields are applied (see also Figured]) and h > is measured in units of J. 

Studies performed by using this confinement geometry show that for certain values of the temperature and the 
competing magnetic fields h, applied to opposite corners, a corner-filling transition can be observed (see Figured]). 
The study of this filling transition under equilibrium conditions has recently attracted growing attention [138-167]. 
Also, the filling transition upon the irreversible growth of a magnetic system has very recently been studied [168]. In 
both cases, the occurrence of an interface between magnetic domains of different orientation is due to the presence of 
competing fields. The localization-delocalization transition of the interface in a finite system yields to a true second- 
order corner-filling transition in the thermodynamic limit (L — * oo). The analytical expression of the equilibrium 
phase diagram was early conjectured by Parry et al. [150] and more recently proved rigorously by Abraham and 
Maciolek [153], yielding 

cosh(2ft/3) = cosh(2A') - e~ 2K sinh 2 (2if ). (9) 

The critical curve obtained by using equation J9]) is shown in Figure 02 It is found that for a given surface 
magnetic field, the filling transition takes place at a temperature lower that of the wetting transition, except of 
course for h = 0, where both curves converge to the Onsager critical temperature of the Ising model (7c). 

In all cases briefly discussed here, one has that the interplay between confinement, boundary conditions, surface 
fields and temperature leads to the occurrence of fluctuating interfaces, so that the propagation of damage in such 
systems is expected to strongly depend on the properties of interfaces. It should also be noticed that due to the 
equivalence between the Ising model and both a lattice gas and a binary alloy, all the physical situations discussed 
above have a wider field of application, e.g. for the study of simple fluids and in the field of condensed matter 
physics. 

5.2 Damage Spreading. 

5.2.1 Strip Geometries. 

Let us start our discussion of DS in confined samples with the case of the strip geometry (LxM, with L <C M) [54,55]. 
Close to criticality, this Ising magnet exhibits an interesting boundary effect, as already discussed in relation to 
Figures CQ (a) and CO(b). When an initial perturbation is introduced in these configurations, at the critical point of 
the DS transition, one observes a monotonic growth of the damage according to a power-law behavior (see Figure 
[6]) given by 

D(t)<xfi, (10) 

where r] is the dynamic critical exponent. Results obtained from the study of DS in these systems are consistent with 
the fact that the presence of interfaces between magnetic domains enhances the propagation of the perturbation, 
as judged by the values of the dynamic critical exponents: r\ OBC — 0.471(5) and r/ PBC = 0.392(5) [54]. Based 
on the fact that rj OBC > rj PBC , one concludes that, for a given time, the damaged area of the sample is always 
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bigger in samples with OBC as compared with samples having PBC. This result can be explained in terms of the 
fluctuations in the orientation of the spins in the region near the interfaces. In fact, around these regions, which 
are present in samples with OBC (see figure Q] (b)), one has the largest fluctuations that enhance the propagation 
of the perturbation. In contrast, inside of the magnetic monodomains, which characterize the samples obtained by 
using PBC (see figure Q] (a)), the propagation of the damage slows down. 

A similar behavior is observed for the propagation of damage in the Abraham's model (see Figure [7]) [56]. In 
this case, for certain values of the temperature T and the surface magnetic fields applied to the upper and lower 
walls of the lattice (h), a critical localization-delocalization transition of the interface between magnetic domains 
occurs. It is found [56] that the dynamic critical exponent for DS is bigger than in the previous cases, yielding 
r] WT = 0.91(1) (see Figure[7]). This result reflects the fact that the interface between magnetic domains is parallel 
to the propagation of the perturbation, and for this reason one has r] WT > rj OBC > rj PBC . Moreover, after proper 
extrapolation to the thermodynamic limit, the phase diagram of the DS transition can be drawn, as shown in Figure 
[3l It is observed that the DS transition occurs within the non-wet phase of the wetting phase diagram (see Figure 
[3| and consequently it does not coincide with the wetting transition. 

Other interesting result related to the propagation of damage along the interface between magnetic domains can 
be visualized by the study of the damage profiles along the L-direction, defined as: 

M 

PlU, t) = — J2 T ) - S U^ T % > ( n ) 

t=i 

where Sfj(t,T) and S B j(t,T) are the reference and damaged configurations at site I of coordinates {i,j}, as early 
defined in the context of equation (1), respectively. Thus, the damage profile represents the average damage of 
the i — th (i = 1, L) row of the system, that runs parallel to the surfaces where the fields are applied, i.e. the 
M -direction. 

It is found [55] that for zero fields and assuming both PBC and OBC (Figure[8](a)-(b), respectively), the damage 
profiles are essentially flat. This can be explained in terms of the interfaces between different magnetic domains 
which are essentially perpendicular to the film and the overall effect on the profiles is homogeneous. On the other 
hand, in presence of surface fields the situation is qualitatively different (see Figure [51[c)). The damage profiles show 
an important curvature related to the existence of an interface in the M -direction, as shown in Figure EJc) . The 
effect of the walls is to slow down the propagation of the damage, while the interface between magnetic domains 
running along the film enhance the spreading of the damage. 

5.2.2 Corner Geometry. 

DS in the corner geometry, with competing (short-range) magnetic fields acting on the surfaces, shows [59] a behavior 
qualitatively similar to that the observed in the case of the strip geometry (see Figure [9]). However, in contrast to 
the previous case, three different regimes for DS were found, as is shown in Figure [TD1 For short times, one observes 
the healing of the damage created in the bulk of the domains. However, at the critical DS point, a small cluster 
of damaged sites survives close to the interface. During the second time regime, the damage propagates along the 
interface according to a power-law behavior of the form D(t) oc t v , with r\* — 0.89(1). Finally, due to the constraint 
imposed by the corners where magnetic fields of opposite direction meet, the damage no longer propagates along 
the surface but starts to spread slowly into the bulk of the domains. Within this late regime one has D(t) oc V , 
where i]** = 0.40(2) is the exponent describing the spreading of the damage in the bulk [59]. 

As in the case of strip geometry, it is useful to gain further insight into the spatiotemporal propagation of the 
damage. For this purpose, the probability distribution of the distance from the damage zone to the corner (P(Iq)) 
can be measured. The distribution was evaluated along the diagonal of the sample (y-direction in Figure [4]) at 
x = L/2. Figure ITTI shows a summary of the results [59]. 

In these cases, for h <C hf (oc) (e.g. h = 0.20) the distribution is almost flat with two small peaks close to the 
corners. Approaching the transition by increasing the field these peaks develop and become slightly shifted toward 
the center of the sample (e.g. h — 0.21 and h = 0.22 in Figure fTTj) . This double-peaked structure indicates that 
the damage remains bound to each corner with the same probability as expected for the case of the non-wet phase. 
On the other hand, for h ~ hf (oo) (e.g. h — 0.23,0.24 in Figure fTTj) the distribution becomes a Gaussian centered 
along the middle of the sample. The Gaussian structure of P(Iq) remains even for h ^> hf (oo) (e.g. h = 0.30 in 
Figure [TT|l . This results are in qualitative agreement with the corresponding to the damage profile showed in Figure 
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[8] for the strip geometry; and it shows that the damage is located in the neighborhood of the interface between 
magnetic domains. 

5.2.3 Overview of the Behaviour of DS in Different Confinement Geometries. 

Table I summarizes the exponents obtained for different confinement geometries and boundary conditions. In view 
of the exponents obtained, one concludes that the presence of interfaces between magnetic domains of different 
orientation causes the enhancement of the propagation as compared with samples where such interfaces are not well 
defined. Furthermore, these results clearly show that the propagation of the damage is anisotropic: it propagates 
better along an interface than across it. Moreover, the obtained results suggest that the damage propagates along 
the magnetic interfaces with a critical exponent rj 1 = 0.90(2) ~ r/ WT ~ rj* and it spreads into the magnetic domains 
with an exponent r] DOM = 0.40(2) ~ ^ BC ~ rj** . 

Taking into account this point of view, the value of the exponent r\ OBC = 0.471(5) can be explained in terms of 
the equilibrium configurations of the system. In fact, the sequence of domains with spins of opposite sign implies 
the inhomogeneity in the propagation of damage. When the perturbation arrives at an interface, it will propagate 
faster than in the bulk. For this reason, the exponent r\ OBC can be thought as a result of the interplay between the 
propagation along the interfaces and the propagation inside the magnetic domains. Regrettably, for this example 
it seems to be impossible to separate both effects in order to write the exponent in terms of and r] DOM . 

Figure [12] shows the behavior of D(t) for the different confinement geometries studied: the strip geometry with 
PBC's and OBC's at criticality and without applied magnetic fields [54,55], the Abraham's Model [56] and the 
corner geometry [59]. In all these cases and because of the dependence on the initial conditions, the first 100 mcs 
were disregarded in order to evaluate the exponents. 



Model 


Critical Exponent rj 


Ising strip geometry with PBC [54] 


0.392(5) 


Ising strip geometry with OBC [54] 


0.471(5) 


Abraham's Model [56] 


0.91(1) 


Corner Geometry (Second Regime) [59] 


0.89(1) 


Corner Geometry (Third Regime) [59] 


0.40(2) 



Table I. List of the critical exponent of Damage Spreading (rj) obtained for the different models, as indicated 
in the first column. 



6 Conclusions. 

A brief overview of the Damage Spreading method applied to different models is presented and the critical non- 
equilibrium DS transition is analyzed. Within this context, DS in confined geometries is discussed in detail. A 
power-law behavior is found for DS as a function of time. The dynamic critical exponent r\ depends on the presence 
of interfaces between magnetic domains generated by the existence of surface magnetic fields applied to the walls 
of the lattice. The results reported suggest that the damage propagates along the interfaces between domains with 
opposite orientations of the magnetization with an exponent rj 1 ~ 0.90 and it spreads into the magnetic domain 
with an exponent rj DOM ~ 0.40. In view of these results, we conclude that the presence of interfaces enhances the 
propagation of a perturbation. 

The critical propagation-nonpropagation transition of the damage was studied too. The results discussed indicate 
that, in some cases as in the Abraham's model and a pure Ising system in the absence of magnetic fields, there 
is clear evidence indicating that the DS transition and the critical transition of the corresponding system occur 
at different critical points. These cases correspond to the wetting and the ferromagnet-paramagnet transitions, 
respectively. In the case of the corner geometry, both transitions —damage spreading and corner-filling— coincide 
within error bars. 

In the cases of confined geometries, the critical exponents found suggest that the DS transition does not belong 
to the DP universality class. This finding seems to be related to the enhancement of the propagation of the damage 
caused by the magnetic interfaces appearing in confined samples. 

The DS method is a suitable tool for the study and understanding of the propagation of perturbations in physical 
systems. Due to the concept of universality that allows one to describe second-order phase transitions and classify 
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them into universality classes, characterized by the same set of critical exponents, it is possible to investigate simple 
models representing more complex physical situations, as e.g. the Ising model as archetype of various systems 
such as ferromagnets, fluids, lattice gases, binary alloys, etc. Therefore, the discussed role of the interfaces in the 
propagation of perturbations may be a quite general phenomenon of relevance in material science in general, and 
particularly, due to the additional influence of confinement, in the field of micro- and nanotechnology. 
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Figure 1: Snapshot configurations corresponding to the Ising magnet as obtained for L = 12, M = 1200, t = 10000 
mcs, T = 0.98Tb, and using (a) periodic boundary conditions, (b) open boundary conditions. Note that the 
horizontal coordinate has been reduced by a factor of five in comparison with the vertical coordinate, for the sake 
of clarity of the picture. Sites taken by down spins are shown in black while up spins are left white [54]. More 
details in the text. 




Figure 2: Snapshot configurations obtained after 10 4 mcs in confined geometries when short-range fields with 
opposite signs are applied to the bottom and top sides of a lattice of size L — 24 and M — 1200. In these cases, the 
snapshots are taken at T = 0.80 Tq and different surface fields: (a) h = 0.4, within the non-wet phase; (b) h = 0.6, 
near the critical wetting curve, and (c) h = 0.8, within the wet phase. Note that the horizontal coordinate has been 
reduced by a factor of five in comparison with the vertical one, for the sake of clarity of the picture. Sites taken by 
down spins are shown in black while up spins are left white [56]. More details in the text. 
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Figure 3: Phase diagram of h versus T/Tc- The dashed line corresponds to the wetting transition (equation ([7])), 
while the full line corresponds to the corner-filling transition (equation ([9])). The circles and squares correspond to 
the results obtained for the DS transition at the strip [56] and the corner [59] geometries, respectively. 
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Figure 4: Corner geometry of size L x L. The signs + and — indicate the surfaces where the competing surface 
magnetic fields are applied. The positive domain is shown in grey and the negative one is left white. In this case, 
the boundary conditions are open for all sides of the sample [59]. 




Figure 5: Snapshot configurations obtained after 10 4 mcs in a corner geometry of size L x L, with L = 128, 
T = 0.80T C and different surface fields: (a) h = 0.20 < h F (L); (b) h = 0.24 - h F {L), and (c) h = 0.28 > h F {L). 
Sites taken by down spins are shown in black while up spins are left white. More details in the text. 
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Figure 6: Log-log plots of D(t) versus t obtained for the Ising magnet in the strip geometry. Results corresponding 
to different temperatures obtained using lattices of size L x M — 12 x 601 and applying open boundary conditions. 
The slope of the dashed line is rf OBC = 0.47 (see equation fT0|) [55]. More details in the text. 



19 



"TTTT 1 1 — I I I I I I I 1 1 1 I I I I I I 1 1 — I I I I I I 



10 



10 



D(t) 



10" 



10" 



• T=0.75 T c 
■ T=0.76 T c 

♦ T=0.78 T c 
▲ T=0.79 T c 
4 T=0.80 T c 
O T=0.82 T c 
n T=0.86 T r 



□ ,0 
n □ / 

□ .0 



1 igisiS 



▲ 



□ 

□ Q0 

X 



□ D p.- 



5~1 



♦ 



10" 



_i i i i i i i i i i i i 



10" 



10' 



1(T 



io J 



l(f 



Figure 7: Log-log plots of D(t) versus t. Results obtained at different temperatures, applying short-range fields of 
magnitude h = 0.5, and using lattices of size L x M = 24 x 1201. The dashed line has slope r) WT = 0.90 [56]. 
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Figure 8: Plots of the damage profile measured along the L-direction and obtained at different times as indicated. 
Results are obtained using lattice of size L = 12, M = 600 in confined geometries, namely (a) T = 0.98Tc and 
PBC, (b) T = 0.98T C and OBC, (c) T = 0.861T C and applying short-range fields = \h L \ = 0.5. [55] 
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Figure 9: Log-log plots of D(t) versus t, obtained for the Ising magnet in the corner geometry and for different 
values of the surface magnetic field h. Results corresponding to T = O.SOTq and obtained using lattices of size 
L = 256. The dashed line has slope rj** = 0.40 [59]. More details in the text. 
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Figure 10: Log-log plot of D{t) versus t, as obtained for the Ising magnet in the corner geometry. Results obtained for 
T = O.SOTc and at the "critical" size-dependent magnetic field hu(L): ho{L = 256) = 0.355, hn(L = 512) = 0.3675, 
h D (L = 1024) = 0.37, and h D (L = 2048) = 0.38. The dashed line has slope rf = 0.89 and the dotted line has slope 
jf* = 0.40 [59]. 
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Figure 11: Plot of probability distribution of the position of the damage P(Iq), obtained for L — 128, T = 0.90Tc, 
and different values of surface magnetic field h, as listed in the figure [59]. 
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Figure 12: Log-log plot of D{t) versus t obtained at the DS critical point for the different models described in this 
work. The dotted line has slope rj DOM — 0.40, the dashed line has slope rj OBC — 0.47, and the solid line has slope 
i] 1 = 0.90. More details in the text. 
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